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SUMMARY 


A theory is presented for calculating stresses in the covers of 
"torsion boxes containing large cutouts. Half the symmetrical uncut por- 
tion of the cover is represented by three stringers that carry all direct 
stresses and the intermediate cover material which is assumed to carry 
only shear stress.' Differential equations of equilibrium are derived 
and solved to obtain stresses along each stringer. Illustrative examples 
are solved and the results are compared with experimental values. For 
one of these examples the results are also compared with a more detailed 
solution made by a numerical method of analysis. The agreement between 
theory and experiment is satisfactory in all cases except those with very 
large cutouts and flexible bulkheads. 


INTRODUCTION 


An approximate method was developed in reference 1 for the analysis 
of stresses in torsion boxes containing large cutouts. This method con- 
sists in first idealizing the structure by concentrating all the direct- 
stress- carrying material of the box into main flange members and assuming 
all the shear to be carried by the webs and cover sheets. Expressions 
are then written for shear and direct stress in each of the members in 
terms of an unknown distribution of torque between the shear webs and 
cover sheets in the cutout bay. Finally, the principle of least work is 
used to find the actual distribution of torque and thus to determine the 
shear and direct stresses. The analysis is adequate for the design of 
all components of such structures except the cover containing the cutout. 

The present paper extends the theory of reference 1 to permit a more 
detailed calculation of the stresses in the covers of such boxes. The 
method is based upon an idealization of the cover similar to that just 
described. except that the full-width portion adjacent to. the cutout is* 
made to contain three stringers on each side of the longitudinal center 
line instead of only one as in reference 1. Differential equations of 
equilibrium are derived and solved to obtain stresses along each stringer. 
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Numerical results for an illustrative example are compared, with experi- 
mental data and also with the results of a solution made by a numerical 
procedure (reference 2). 
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SYMBOLS 

length of gross section at one end of box, inches 

half- width of box, center to center of webs, inches 

(bp + 2b 2 ) 

width of net section, inches 

width of each of the two interior panels in gross section 
of idealized structure, inches 

depth of box, center to center of cover sheets, inches 

half-length of net section 

fraction of total torque carried by shear webs in cutout 
bay 

thickness of cover sheet, inches 

thickness of top cover sheet in net section, inches 

thickness of top or bottom cover sheet in gross section, 
inches 

thickness of lower cover sheet in net section, inches 

thickness of torque-transfer bulkhead, inches 

thickness of shear web in gross section, inches 

thickness of shear web in net section, inches 

displacements in the x-direction at any point along 
stringers 1, 2, and 3, respectively, inches 

displacement of any cross section in y-direction, inches 


x, y 


coordinate axes 
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y , y distance from neutral axis of net section to stringers 1 

• and 2 , respectively, inches 

A-^, Ag, A^ areas of stringers 1, 2, and 3, respectively, of idealized 
cover, square inches 

A^, A^. . .A 7 areas of flanges in simplified box, square inches (see 

fig. 7 'for location of each area) 


^i* Cg. . .c 6 

D 

E 

G 

I 

K, K* 


constants of integration 
differential operator 

Young's modulus, pounds per square inch 

shear modulus, pounds per square inch 

moment of inertia of net section, inches 

dimensionless constants leading to the solution for the 
parameter k 


T 

PiJ 




torque, inch-pounds 

factors in assumed solutions of differential equations 
(i designates the stringer and j designates the 
characteristic value in the solution of the differential 
equations of equilibrium) 

characteristic values in the solution of the differential 
equations of equilibrium 


a l’ a 2’ a 3 direct stress in gross • section in stringers 1, 2, and 3^ 
respectively, pounds per square inch 

crj , <Tg direct stress in net section in stringers 1 and 2, respec- 

^ n tively, pounds per square inch 


T l> t 2> t 3 


shear stress in gross section in panels 1, 2, and 3, 
respectively, pounds per square inch 


'n 


shear stress in net section, pounds per square inch 


METHOD OF ANALYSIS 


, The problem considered is the analysis of the cover of a torsion box 
of rectangular cross section with a centrally located cutout as indicated 
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in figure 1. As a first step in the analysis, the cover containing the 
cutout is separated from the rest of the box as a "free body." In place 
of the portion of the box that has been removed, forces must be added to 
the cover that are equal to those exerted by the box on the cover -when 
the two parts were fastened together. These forces may be calculated by 
the method of analysis given in reference 1. Such a procedure leaves a 
stringer-stiffened panel, subjected to forces as indicated in figure 2, 
to be analyzed. Along the sides of the panel are forces representing the 
loading transferred to the corner angle by the shear webs. The lateral 
forces along the ends come from the end bulkheads, and the additional 
lateral forces at the junction of the full-width and cutout sections of 
the cover are introduced by the torque transfer bulkheads that bound the 
cutout bay. If the panel is assumed to be doubly symmetrical, only one 
quarter of it, as indicated by the shaded area in figure 2, need be con- 
sidered. As before, the portion of the cover that is cut away must be 
replaced by the forces it exerts. The rest of the cover is an L- shaped, . 
stringer-stiffened panel subjected to forces as shown in figure 3. 

To simplify the analysis further, the panel is assumed to be repre- 
sented with sufficient accuracy by the idealized structure of figure k, 
which consists of three stringers assumed to carry all the tensile and 
compressive stress and three panels of cover sheet assumed to carry all 
the shear. This simplification follows a line of reasoning similar to 
that used in reference 3 . The three stringers are located as follows: 
Stringer 1 is placed along the line of intersection of the shear web and 
cover sheet. Stringer 2 is placed along the edge of the cutout. Stringer 3 
is placed midway between stringer 2 and the longitudinal center line of the 
cover. Between the stringers is a flat sheet of the same thickness as the 
cover sheet of the actual box but which differs from it by the assumption 
that it can carry only shearing stress. All the direct-stress-carrying 
material in the cover is concentrated in the three stringers. Distribu- 
tion of this 'material to the different stringers is assumed to be as 
follows: All the corner angle is included in stringer 1. All' the coaming 

stringer (the stringer bordering the cutout) is included in stringer 2. 

The rest of the material, including cover sheet and stiffeners, which lies 
between the locations of stringers 1 and 2 is divided equally between them. 
All the rest of the cover is included in stringer 3* 

. The analysis of the idealized structure may be subdivided further 
by cutting the frame along the dashed line in figure 3 . The narrower 
portion that borders the cutout is referred to as the net section, and 
the full-width portion is called the gross section. The gross section 
and the external forces acting on it are shown in figure 5 and similar 
information for the net section is given in figure 6. . 

The factor k is taken from reference 1. The various dimensions 
entering into its calculation are shown in figure 7 and the method for 
computing it is given in appendix A. 
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Stresses in the net section may he calculated with sufficient 
accuracy by the ordinary theory of cantilever beams. Thus the stress 
in the edge stringer is equal to 


' “In * - *> - 

and the stress in the coaming stringer is given by 

T y 2 , 

(Jo = — — (l - k) x 
■• n 2cl 


( 1 ) 


( 2 ) 


The shear stress is constant throughout the net section and is equal to 


T. 

2b^ct 


(1 - k) 


(3) 


Stringer stresses in the gross section may be computed from the 
equations 


where 




> ( 6 ) 


2, 3) 


(j = 1 , 
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where 


h - * |‘i[ 2 - »2(^)\j 2 ] + Asf - b ;=(^>j 2 ] f - T • b ^)^ 
A 3 [I- h. 

(use whichever sign is necessary to make p 2 j positive) 


Gt 


1 > 2 EA 1 A 2 A 3 hh 1 h 2 


^k^y^bb^ + i'AjA^b^b^ + b^ + Ar>A 3 b 2 i 


|Ai 2 A 2 2b2 bi 2 -+ i A 1 2 A 3 2 b 2 (b 1 + b^) 2 + A^A^b^ - 
A 1 A 2 A 3 bb 2 (A 1 bb 1 + 2A 2 b ] b 2 - Agb^ 


r (7) 


X 3 2 = 0 


(X^ is associated with the positive root). The coordinate x is measured 
from the junction of the net and gross section toward the tip. Shear 
stresses are computed from the equations 


ri - 55t E - ^ ^ + T 

r 2 = T 1 + T 5? 


A 1 d£i 
dx 




Ao dCq 

T 3 = T 2 + t dx 


Details of the derivation of the equations are given in appendix B. 


(8) 
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NUMERICAL EXAMPLE 


In order to demonstrate the computational procedure, the method is 
applied to the torsion hox designated as case 1 in reference 1. After 
isolating the cover from the rest of the hox and idealizing it in the 
manner outlined previously, the following data are obtained: 

Ap = 1.936 square inches A2 = 0.575 square inches 

A ^ = 1.306 square inches 

h = 25.70 inches bp = 11.12 inches 

h 2 = 7*29 inches . 

a = 36.75 inches c = 10 inches 

d = 24.5 inches t = O.O63 inches 

I = 36.75 inches 1 * T = 99 >750 pound-inches 

y-j_ = 2.30 inches = 8.82 inches 

The shear and Young's moduli are taken as 

G = 4,000,000 pounds per square inch 
E = 10,600,000 pounds per square inch 
The value of k, obtained from reference’ 1, is 

k =0.546 

Stresses in the net s'ection are then computed from equations (l), 
(2 ), and (3) as 

Op^ = 87x 
a 2 n = 543x 
T n = 3,230 
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The first step in computing the direct stresses in the gross section 
is to solve equations (.7) for Xg, and X^. The results are 



x x = 0.1088 

X2 

= 0.0564 

x 3 = 0 

Equations (6) 

are solved for p^j 

and yield 

the following values: 

3 ll 

= - 0.1124 

P21 = 

1.1413 

p 31 = - 0.4165 

Pl2 

= - 0.2340 

p 22 = 

0.5412 

P32 = 0.7454 

P13 

= 0.6702 

P23 = 

0. 3802 

P33 = 0.1901 

Values of the 
as follows: 

constants 

c 2> 

and are Calculated from equation (5) 


Ci = .5963 

C 2 

= -3260 

c 3 = -293 


Substitution of the values for and the integration constants into 

equation ( 4 ) yields 


a 1 = 670e -0,1088x + 763e"°*° 56 ^ x - 196 + 5-3x 
a 2 = -6805e“°- lo88x - 1764e"°‘° 56l *' x - 111 + 3. Ox 
a 3 = 2483e"°* 1088x - 2430e"°*° 56 ^ x - 56 + 1.5x ■ 

Shear stresses are computed from equations (8) and the results are 
t 1 = 1610 - 224le"°' 1088x - 1322e -0,056 ^ x 
t 2 = 1637 + 4517e" 0-1088x - 4l4e“°*°^ x 

t 

t 3 = 1669 - 1084e" 0,1088x + 2427e -0 * ®^64x 



NACA TN 2290 


9 


Figures 8 and 9 show the shear and stringer stresses, respectively, 
computed from these equations. Similar calculations have been made for 
the other cases listed in reference 1, and the results are shown in fig- 
ures 10 to 15 . For convenience of reference, the necessary data for 
each case are given in table 1. Cases 1, 2, and 3 differ only in the 
width of the cutout. Case 3 bas practically a full- width cutout. Case k 
has been omitted here, since it differs from ease 3 only by the removal 
of intermediate bulkheads, which. does not affect the computed values. 

Case 5 bas a slightly longer cutout than case 3 > but of the same width. 
Principally, however, it differs from case 3 in that the torque- transfer 
bulkheads are made very much more rigid. 


RESULTS AND DISCUSSION 


The preceding analysis is based upon an idealization of the structure, 
which permits the calculation of direct stresses along only three longi- 
tudinal lines. Between stringers, the variation of direct stress is 
assumed to be linear. The cover is divided into three panels, and the 
shear is assumed to be constant in each panel at any cross section. In 
addition, the forces which are assumed to be acting on the cover also 
are calculated by an approximate method based upon an idealized structure. 
These considerations make it clear that only approximate results may be 
expected and that the stress distribution can be indicated in only a very 
general way. Upon such a basis, the correlation between ' computed and 
measured stresses may be considered good. The theory of reference 1 is 
capable of predicting only a constant shear stress throughout the gross 
section of the cover, whereas the experimental results indicate that the 
true distribution deviates greatly from a uniform stress, especially in 
the vicinity of the cutout. The method outlined in the present paper 
permits the shear stresses to be divided over three regions in each half 
of the cover, represented by the three panels between the stringers, and, 
in addition, the stresses may be varied from one cross section to another. 
The true state of stress can therefore be approximated much more closely. 

In a similar way, the additional stringers make possible a more detailed 
distribution of the direct stress. The theory appears to indicate peak 
shear values a little below those actually existing, as would be expected, 
since the constant shear represents average stress over the panel. 

Stringer stresses seem to be fairly reliable in the case of cutouts of 
moderate width or those that are bounded by very rigid bulkheads, but 
for cutouts of almost full width combined with flexible bulkheads, some 
sizable discrepancies occur. For instance, cases 1, 2, and 3 all have 
quite flexible bulkheads. Case 1 represents the smallest cutout for 
which calculations were made, and the experimental data follow the com- 
puted curves quite closely, as indicated by figure 9* Case 2 represents 
a somewhat wider cutout, and the stringer stresses begin to show some 
deviation near the cutout as shown at station 27«5 in figure 11. Case 3 
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represents almost a full-width cutout, still with a flexible bulkhead, 
and the stringer stresses show a marked discrepancy near the cutout as 
at station 27.5 in figure 13. Case 5 shows practically the same cutout 
as case 3 but with a greatly strengthened torque- transfer bulkhead. 

Computed and measured stringer stresses agree much more closely as 
indicated by figure 15. Case 3 is somewhat extreme and for such structures 
it may be necessary to develop a more comprehensive theory which takes 
account of additional factors such as distortion in the chordwise direction. 

Case 1 was computed by a numerical method discussed in appendix C. 
Results are shown by the dashed lines in figures 8 and 9. The method 
is somewhat more accurate than the one already described but was judged 
not to offer sufficient advantages to justify the added labor that it 
requires. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., October 27, 1950 
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APPENDIX A 
COMPUTATION OF k 


The external forces that are shown in figures 3 and. 5 have been 
derived from reference 1 and are written as functions of k. The 
factor k represents the fraction of the total torque carried by the 
shear webs in the cutout bay. No derivation of k will be given here, 
but for convenience, the equations necessary for computing it are as 
follows: 

10 

X K ’ 

■ k.-i 

10 

X> 

1 • 

where 
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k 6* 



_ 1*G b 2 d 2 
f - 3E o 

Vi 


k 7 ' 



^ kG d 2 / 

K8 sx;r 



K, 


1 _ 


■8 = 


b + 2bn 



16 G 

3E a 6 



v _ 32G ad 

10 " ^“.a 7 


k io' “ 2 K i° 


Subscripts for the various A terms are as indicated in figure J. 
Locations of the stringers having areas Aj^ to Aj inclusive also 
are shown in the same figure. The amount of area to include in the 
idealized stringers is a matter of judgment to some extent. For the 
structures considered here, the choice is made relatively easy because 
of the large corner angles. Areas A^ to Aj inclusive may be taken 

as only the comer angle, and all other areas neglected. In other cases, 
however, ' consideration must be given to the contributions of the shear 
webs, cover sheets, and stringers. The area of the coaming stringer and 
a suitable addition for the influence of the rest of the cover 6f the 
net section is included in A^. 
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APPENDIX B 

DERIVATION OF STRESS EQUATIONS FOR GROSS SECTION 


Figure 5 shows the idealized gross section and the external loading 
system to which it is subjected. Expressions for the loads are taken 
from reference 1, except for the two loads applied to stringers 1 and 2 
which are transmitted by the stringers of the net section. These loads 
are derived from static considerations of the net section. With the 
notation 


U 1 

u 2 

u 3 


displacement of any point 
displacement of any point 
displacement of any point 
displacement of any cross 


on stringer 1 in the x-direction 
on stringer 2 in the x-direction 
on stringer 3 in the x-direction 
section in the y-direction 


and the positive shear defined as that which distorts an element of area 
so that the length of the first- quadrant diagonal is increased, the 
following relationships may be written: 

dui cr-^ 
dx E 


du 2 ag 

d 7“ E 

du^ 03 
dx~ “ IT 



(Bl) 







I 1 * . 
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Summing forces in the x-direction along each stringer in turn yields 






(B3b) 


(B3c) 
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In order to obtain an expression for dv/dx, use is made of the fact 
that the sum of the internal shear forces must equal the external load, 
or 


j£[l + |(2k - 1)] = t + T 2 b 2 + T 3 b 2 ) 


(B4)- 


Substituting the shear forces given in equations (Bl) into equation (B4) 
and simplifying yields 


or 


dv 

dx 



— F 

l*Gtbc _ 


+ -(2k - 



> 


d^v _ °1 

dx^ Eh 


(B5) 


Upon substitution of the equation for d^v/dx^ into equations (B2) and 
then (B2) into (B3), the following differential equations are obtained: 
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Introducing the notation = D 2 a and collecting terms allows 

dx 2 

the differential equations (B6) to be written* in matrix form as 


A X ED 2 - Gt(— - 

W b > 

Gt 

bl 

0 


a l 

Gt 

*1 

AoED 2 - Gt (— + — ) 

^ V b l b 2/ 

Gt 

b 2 


0 2 

0 

Gt 

AoED 2 - — 

3 b 2 


a 

b 2 


3 


= 0 (B7) 


To solve the differential equations, assume solutions of the form 

A.x 


on on 

^ = Cp x e a 2 = CP 2 e 


a ^ = Cp^e 


where C = Constant of integration.' Then 


D 2 a 1 = CPjX^* 


D c a 2 = CP 0 X e 


2 Xx 


and equations (B7) become 


A X EX^ - Gt 


ft - e 


D 2 ^ = Cp^X^*’* 


Gt 

b. 


Gt 

*1 


AgEX 2 - Gt 


(£- + i 

V 5 ! *2. 


Gt 


Gt 


A 0 EX 2 - — 


b 2 


= 0 (b8) 
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from which, "by setting the determinant of the coefficients equal to zero, 
2 7 
three values of X are found 


%2 


Gt 




^A^Agbb^ + i AjA^b + ^2^ 3^2 " 


Ai^A^^b^^i 4. + A2?A2^b2^ ™ 

, gTl 1 / 2 ! 

^1^2^3^2v^l^l ^ 2k^o^2 ~ ^3^2 /I 


> (B9) 


2 

X 3 =0 


(X-^ is associated with the plus root). There is a separate set of values 
of 3 corresponding to each different value of X^. These values of 3 
are denoted hy 3ii> P-^j and Pj^j the second subscript corresponding 

to the subscript of the Xj 2 with which it is associated. The values of 
the various 3^ may be determined in relation to each other from two of 
the equations (B8). Thus 


3 U 


ii 


1 - - b. 


EA 


— “ ^1 2 
Gt J 


3 


3j 


2j 


EA, 2 
2 - b^ — X .• 


' 2 gT J 


> (BIO) 


(J = 1, 2, 3) 
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The following expressions for may now be written: 


_ ^lPll e ^" X + ^2^12 e + ^ 3^13 + *-'4Pl3 x + Ccj(3-j^e^l x + C6Pp2 e ^ X 

a 2 = C l |3 21 e ^ lX + C 2^22 e + C 3 P 23 + C M 3 23 X + + c 6^22 e ^ 2X L(B11 ) 


a 3 C l p 31 e 1 + C 2 P 32 e 2 + c 3 p 33 + C 4 P 33 X + C 5 P 31 e 1 + C 6 P 32 e 2 


This result constitutes the formal solution of the differential equations. 

The integration constants C-^ Cg . . . Cg are as yet arbitrary 

and are used to make the stresses satisfy the boundary conditions. The 
boundary condition at each end of each stringer is substituted into the 
stress equations (Bll) and the resulting six equations are solved simul- 
taneously for C-l, C 2 . . . C 5 . The integration constants can be 
evaluated much more easily, however, by taking advantage of an orthogo- 
nality relation which may be shown to exist among the values of 

Before proceeding with this step, it may be noted that, if the box 
is infinitely long, C<p and Cg must vanish in order to obtain finite 
stress at the tip. Calculations for a number of boxes of practical pro- 
portions indicate that only a negligible error is introducted by dis- 
regarding these two terms. Thus, for practical purposes, the stress 
equations (Bll) may be reduced to 


a ± - ^]_Pii e ^ + ^2^i2 e ^ + ^3^13 + C 4^13 x 

(B12a) 

~A,,x -A. x 


°2 = ^1^21 e ^ + ^2^22 e ^ + ^3^23 C 4^23 x 

(B12b) 

a 3 = C-^P^i e ^ + ^2^32 e ^ + C3P33 + C^33X, 

(B12c) 
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At x = 0, the boundary conditions are 


2Apbpc 


k - k 


a 9 = (k - 1) 

2 2A 2 bpC ' 


°3 ' 0 


Substitution of these values into equations (B12) gives 


CpPpp + ^2^12 + ^ 3^13 = 


2A^b-^c 


E - 4 + t)] 


C 1 P 21 + C 2 P 22 + C 3 P 23 2A 2 b lC ^ k ~ 


C1P31 + C2P32 + C 3^33 ” 0 


It may be shovn by manipulation of equations (B8) that 


^■ 1 ^ 11^12 + ^ 2 ^ 21^22 + ^ 3 ^ 31^32 = ® 


A 1 P 12 P 13 + A 2 P 22 P 23 + A 3 P 32 P 33 ° 


A l p ll p 13 + A 2 P 21 P 23 + A 3 |3 31 |3 33 " 0 


For instance, if equations (B8) for Xj' = are multiplied by 0 12 , 
equations (B8) for kj = A, 2 ' are multiplied by Ppp, the corresponding 
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equations of each set are subtracted from .each other, and finally the 
resulting three equations are added, then the first of equations (Bl4) 
is obtained. A similar procedure yields the remaining two equations. 


Also, and have been expressed as 

.(equations (BIO)) which is as yet unrestricted, 
be obtained by choosing 80 that 


functions of 
Some simplification' may 


+ ^2^22^ + A3P33 2 = 1 


(B15) 


in a manner analogous to normalization. 


Now, let equation (B13a) be multiplied by A 1 .p 11 , equation (B13b) by 
A 2^21-> and elation (B13c) by A^P^p and let the results be added. If 

the orthoganality relations (Bl4) and equation (B15) are taken into 
account, it is found that 


Cl '^{[- k ( 1 -Tj 11+<k - 1)e 


In a similar manner, it may be determined that 


21 | 


and 


c 2 --SM 

* 2 bpc 


„ Td 
Co - < 

J 2 b x c 


i - k ( 1 + S] 

i - k ( 1 + i§ 


p 12 + <k 


1)0 


221 


P 13 + (k - 1)3 


231 


(B16) 


(B17) 


(B18) 


At x = a the boundary conditions are a p = a 2 = a 3 = With 

these values substituted into equations (B12), equation (B12a) multiplied 
By A l3l3> equation (B12b) by A 2 3 2 ^, and equation (B12c) by A 3333 , 

addition 'of the equations yields 



(B19) 
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Introduction of the values for C 2 , C and C]^ into equations (B12) 

defines final solutions for the stress in each of the three stringers at 
any point x distant from the y-axis. 

Shearing stresses in each of the three panels may be obtained 
directly from equations (B3) as 


■ 1= * C-ata-yl + ii£i 

4bct I— a _J t dx 


t 2 = T i + 


A 2 dcr, 


> 


t dx 


(B20) 


Ao dOq 

t 3 = t 2 + ~T — ~ 
^ ^ * dx 



22 


NACA TN 2290 


APPENDIX C 

DISCUSSION OF NUMERICAL PROCEDURE 


The idealized cover used in the preceding method could he made more 
nearly like the structure from which it was derived by increasing the 
number of stringers. It is clear that such a procedure would greatly 
increase the labor necessary to obtain a solution, especially since the 
equation for X would be of some higher degree. If additional stringers 
must be included in the analysis, it probably is simpler to use a numerical 
procedure such as the one given in reference 2. In this method, the 
external forces acting on the cover are computed in the same manner as 
before, but the idealization of the structure takes a different form. 

The cover is divided into a number of rectangular units or panels as 
indicated in figure l6. Each of these panels is assumed to be bounded 
on the sides parallel to the longitudinal center line of the box by 
stringers of finite extensional stiffness. The panels also are assumed 
to be bounded on the ends by infinitely rigid members which maintain the 
cross section but offer no resistance to motion outside of that plane. 
Between the edge members, it is assumed that a sheet of the same thick- 
ness as the cover exists and that it carries only shearing stress; 

If equilibrium of the elements of the structure is to be satisfied, 
a definite relationship must exist between the displacement of any given 
panel and the displacements of all the panels that form its borders. 

This relationship depends upon the elastic properties of each of the 
panels involved. Details of these relationships and their derivation 
are given in reference 2. When such relationships (hereinafter referred 
to as iteration equations) have been established for all the panels, an 
assumed set of displacements is introduced. The only restriction upon 
these displacements is that they must satisfy the statics of the problem, 
but, of course, the work involved is reduced if the assumed values are 
as close as possible to the actual displacements. A usual assumption is 
that the values will conform to the standard rules of strength of materials. 
If all the assumed displacements were exactly correct, all the iteration 
equations would be satisfied. If the assumed displacements are not exact, 
new values for the displacements are obtained which more nearly satisfy 
the iteration equations. The process may be repeated until sufficiently 
accurate displacements are obtained. 

\ 

When the iteration has been completed, the results are in the form 
of displacements of each comer of each panel into which the cover was 
divided. With this information, it is a simple matter to calculate shear 
or normal stresses along any desired section. The detail with which the 
stress distribution can be obtained depends upon the number of panels 
which are included in the idealized cover. 
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The same problem that was used as an illustrative example of the 
three- stringer method has been solved by the numerical procedure. The 
manner in which it was subdivided into panels is illustrated in figure 16. 
Shear and normal stresses are indicated by the dashed lines in figure's 8 
and 9- 
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